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ABSTRACT 

We report the discovery of two new transiting extrasolar planets. HAT-P-18b orbits the V =12.759 
K2 dwarf star GSC 2594-00646, with a period P = 5.508023 ± 0.000006 d, transit epoch T c = 
2454715.02174 ± 0.00020 (BJD), and transit duration 0.1131 ± 0.0009d. The host star has a 
mass of 0.77 ± 0.03 M Q , radius of 0.75 ± 0.04 i? Q , effective temperature 4803 ± 80 K, and metal- 
licity [Fe/H] = +0.10 ± 0.08. The planetary companion has a mass of 0.197 ± 0.013Mj, and 
radius of 0.995 ± 0.052 Rj yielding a mean density of 0.25 ± 0.04 gem" 3 . HAT-P-19b orbits the 
V =12.901 Kl dwarf star GSC 2283-00589, with a period P = 4.008778 ± 0.000006 d, transit epoch 
T c = 2455091.53417 ± 0.00034 (BJD), and transit duration 0.1182 ± 0.0014d. The host star has a 
mass of 0.84 ± 0.04M Q , radius of 0.82 ± 0.05 i? , effective temperature 4990 ± 130K, and metallic- 
ity [Fe/H] = +0.23 ± 0.08. The planetary companion has a mass of 0.292 ± 0.018 Mj, and radius 
of 1.132 ± 0.072 Rj yielding a mean density of 0.25 ± 0.04gcm -3 . The radial velocity residuals for 
HAT-P-19 exhibit a linear trend in time, which indicates the presence of a third body in the system. 
Comparing these observations with theoretical models, we find that HAT-P-18b and HAT-P-19b are 
each consistent with a hydrogen-helium dominated gas giant planet with negligible core mass. HAT-P- 
18b and HAT-P-19b join HAT-P-12b and WASP-21b in an emerging group of low-density Saturn-mass 
planets, with negligible inferred core masses. However, unlike HAT-P-12b and WASP-21b, both HAT- 
P-18b and HAT-P-19b orbit stars with super-solar metallicity. This calls into question the heretofore 
suggestive correlation between the inferred core mass and host star metallicity for Saturn-mass plan- 
ets. 

Subject headings: planetary systems — stars: individual ( HAT-P-18, GSC 2594-00646, HAT-P-19, 
GSC 2283-00589 ) techniques: spectroscopic, photometric 



1. INTRODUCTION 

Extrasolar planets which transit their host stars (Tran- 
siting Extrasolar Planets, or TEPs) provide a unique 
opportunity to determine the bulk physical properties 
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(mass, radius and average d ensity) of planetary bodies 
outside the Solar System (e.g. lCharbonneaull2009D . From 
the more than 90 such planets that have been announced 
to datJ^I. it has become apparent that gas giant plan- 
ets more massive than 0.4Mj exhi bit a wide range of 
radii (from 0.885i?j for C oRoT-13b.lCabrera et al.| [2010l 
to 1.79i?j for WASP-12b. lHebb et al.ll2009D . Below this 
mass fewer planets are known; however the seven known 
TEPs with masses similar to Saturn (0 15M T < M < 
0.4Mj; the mass of Saturn is 0.299 Mj, iStandishl [l995l) 
also appear to have diverse bulk properties. Two of 
these TEPs have densities much less than that of Saturn 
(HAT-P-12b and WASP-21b both have p ~ 0.3 gem" 3 
while Saturn has p ~ 0.7gcm -3 ; lHartman et al.1 [20091 : 
iBouchv et al.l 120101 ) . three have densities that are some- 
what less than that of Saturn (Kepler-9b and Kepler- 
9c have dens ities of p ~ 0.5gcm~ 3 and p ~ 0.4gcm~ 3 
respectively, iHolman et all 120101: and WASP-29b has 
p ~ 0.65 gcm~ 3 . iHellier et al.ll2010l ). and two have den- 
sities that are greater than that of Saturn (HD 149026b 
has p ~ 0.85gcm- 3 . lSato et al.l l 2005t ICarter et aJjl2009t 
and CoRoT-8b has p ~ 1.6 gem" 3 . [Borde et al.l T2010). 
The inferred core masses of these planets also differ dra- 
matically, with the two low-density planets having neg- 
ligible cores of Mq % lOM®, the three intermediate- 
density planets having cores that are perhaps several tens 

12 e.g. see http://exoplanet.eu 
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of Earth masses, and the two high-density planets hav- 
ing cores that represent a substantial fraction of their 
respective masses. The inferred core masses of planets 
in this mass range appear to correlate with the metal- 
licity of the host star. The two low density planets or- 
bit stars with sub-solar metallicity ([Fe/H]= —0.29 for 
HAT-P-12, and [Fe/H]= -0.4 for WASP-21), while the 
five higher density planets orbit stars with super-solar 
metallicity ([Fe/H]= 0.3 for CoRoT-8, [Fe/H] = 0.36 for 
HD 149026, [Fe/H]= 0.11 for WASP-29, and [Fe/H] = 
0.12 for Keplcr-9). This has been taken as suggestive evi- 
dence for the core-accretion scenario for planet formation 
(lAlibert et all 120051; iGuillot et all 12009: lHartman et al.l 
120091: IBouchv et alJl20ldT ~ 

In this work we present the discovery of two new low- 
density planets with masses comparable to that of Sat- 
urn. The new planets HAT-P-18b and HAT-P-19b have 
masses that are very similar to HAT-P-12b and WASP- 
21b respectively, and have densities that are slightly less 
than each of these planets. However, both new planets 
orbit stars with super-solar metallicity, casting doubt on 
the correlation between planetary core mass and stellar 
metallicity for Saturn-mass planets. 

The planets presented in this paper were discovered 
by the H ungarian-made Au tomated Telescope Network 
(HATNet; Bak os et al.ll2004l) survey, which has been one 
of the main contributors to the discovery of TEPs. In 
operation since 2003, it has now covered approximately 
14% of the sky, searching for TEPs around bright stars 
(8 < I < 14). HATNet operates six wide-field instru- 
ments: four at the Fred Lawrence Whipple Observatory 
(FLWO) in Arizona, and two on the roof of the hangar 
servicing the Smithsonian Astrophysical Observatory's 
Submillimeter Array, in Hawaii. Since 2006, HATNet 
has found seventeen TEPs. In this work we report our 
eighteenth and nineteenth discoveries, around the rela- 
tively bright stars also known as GSC 2594-00646, and 
GSC 2283-00589. 

The layout of the paper is as follows. In Section [2] we 
report the detections of the photometric signals and the 
follow-up spectroscopic and photometric observations for 
each of the planets. In Section [3] we describe the anal- 
ysis of the data, beginning with the determination of 
the stellar parameters, continuing with a discussion of 
the methods used to rule out nonplanetary, false posi- 
tive scenarios which could mimic the photometric and 
spectroscopic observations, and finishing with a descrip- 
tion of our global modeling of the photometry and radial 
velocities (RVs). Our findings are discussed in Sectional 

2. OBSERVATIONS 

2.1. Photometric detection 

Table Q] summarizes the HATNet discovery observa- 
tions of each new planetary system. The calibration 
of the HATNet frames was carried out using standard 
photometric procedures. The calibrated images were 
then subjected to star detection and astromctric deter- 
mination, as described in lPal fe B akos (2006). Aperture 
photometry was performed on each image at the stel- 
lar centroids derived from the Two Micron All Sky Sur- 
vey (2MASS: ISkrutskie et al1l2006ft catalog and the indi- 
vidual astrometric solutions. The resulting light curves 
were decorrelated (cleaned of trends) using the External 



Parameter Decorrelation (EPD; see iBakos et al.l I2010T) 
technique in "constan t" mode and t he Tr end Filtering 
Algorithm (TFA; see IKovacs et all 120051 ). The light 
curves were searched for periodic bo x-shaped signals us- 
ing the Box Least-Squares (BLS; see IKovacs et al.ll2002| ) 
method. We detected significant signals in the light 
curves of the stars summarized below: 

• HAT-P-18 - GSC 2594-00646 (also known as 
2MASS 17052315+3300450; a = 17 h 05 m 23.28s, 
5 = + 33°00'45.0"; J2000; V = 12.759 lDroege et all 
2006). A signal was detected for this star with an 
apparent depth of ~18.5mmag, and a period of 
P =5.5080 days (see Figure!]}. The drop in bright- 
ness had a hrst-to-last-contact duration, relative to 
the total period, of q =0.0205 ±0.0002, correspond- 
ing to a total duration of Pq =2.716 ± 0.021 hr. 

• HAT-P-19 - GSC 2283-00589 (also known as 
2MASS 00380401+3442416; a = 00 h 38 m 04.02s, 
5 = + 34°42'41.7"; J2000; V = 12.901 IDroege et all 
2006). A signal was detected for this star with an 
apparent depth of ^22.0 mmag, and a period of 
P =4.0088 days (see Figured]). The drop in bright- 
ness had a first-to-last-contact duration, relative to 
the total period, of q =0.0295 + 0.0004, correspond- 
ing to a total duration of Pq =2.837 ± 0.034 hr. 
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Fig. 1— Unbinned light curve of HAT-P-18 including all 10,000 
instrumental /-band 5.5 minute cadence measurements obtained 
with the HAT-6 and HAT-9 telescopes of HATNet (see Table [TJ, 
and folded with the period P = 5.5080228 days resulting from the 
global fit described in Section [3] The solid line shows a simplified 
transit model fit to the light curve ("Section 13. 311 . The bold points 
in the lower panel show the light curve binned in phase with a 
bin-size of 0.002. 



2.2. Reconnaissance Spectroscopy 

As is routine in the HATNet project, all candidates 
are subjected to careful scrutiny before investing valu- 
able time on large telescopes. This includes spectro- 
scopic observations at relatively modest facilities to es- 
tablish whether the transit-like feature in the light curve 
of a candidate might be due to astrophysical phenomena 
other than a planet transiting a star. Many of such false 
positives are associated with large RV variations in the 
star (tens of kins -1 ) that are easily recognized. We made 
use of three different facilities to conduct these obser- 
vations, including the Harvard-Smithsonian Center for 
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TABLE 1 
Summary of photometric observations 



Instrument/Field 



Datc(s) 



Number of Images Cadence (s) Filter 



HAT-P-18 



HAT-6/G239 


2007 Mar-2007 Jun 


4383 


330 


1 


HAT-9/G239 


2007 Mar-2007 Jun 


5719 


330 


I 


KeplerCam 


2008 Apr 25 


215 


73 


Sloan i 


KeplerCam 


2009 May 10 


185 


133 


Sloan g 


HAT-P-19 










HAT-7/G163 


2007 Sep-2008 Jan 


2324 


330 


R 


HAT-8/G163 


2007 Sep-2008 Jan 


1617 


330 


h> 


HAT-6/G164 


2007 Sep-2008 Feb 


3676 


330 


R 


HAT-9/G164 


2007 Sep-2008 Feb 


2711 


330 


R 


KeplerCam 


2009 Oct 06 


37 


81 


Sloan i 


KeplerCam 


2009 Oct 30 


97 


150 


Sloan i 


KeplerCam 


2009 Nov 27 


76 


81 


Sloan i 


KeplerCam 


2009 Dec 01 


194 


89 


Sloan i 
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Fig. 2.— Unbinned light curve of HAT-P-19 including all 10,000 
instrumental .R-band 5.5 minute cadence measurements obtained 
with the HAT-6, HAT-7, HAT-8 and HAT-9 telescopes of HATNet 
(see Table Q}, and folded with the period P = 4.0087782 days re- 
sulting from the global fit described in Section [3] The soli d lin e 
shows a simplified transit model fit to the light curve (Section 13.31 1. 
The bold points in the lower panel show the light curve binned in 
phase with a bin-size of 0.002. 



Astrophysics (CiA) Digital Speedometer (DS; ILathaml 
1992), and th e Tillinghast Re flector Echelle Spectro- 
graph (TRES; iFuresz. G.II2008I ). both on the 1.5 m Till- 
inghast Reflector at the Whipple Observatory on Mount 
Hopkins, Ariz ona, and the Fibre-fed Ec helle Spectro- 
graph (FIES; iFrandsen fe Lindberdll999f) on the 2.5 m 
Nord ic Optical Telescope (NOT; Diupvik fe Andersen! 
2010) at La Palma, Spain. We used these facilities to 
obtain high-resolution spectra, with typically low signal- 
to-noise (S/N) ratios that are nevertheless sufficient to 
derive RVs with moderate precisions of 0.5-1.0 km s" 1 
for slowly rotating stars. We also use these spectra to es- 
timate the effective temperatures, surface gravities, and 
projected rotational velocities of the stars. With these 
observations we are able to reject many types of false 
positives, such as F dwarfs orbited by M dwarfs, grazing 
eclipsing binaries, or triple or quadruple star systems. 
Additional tests are performed with other spectroscopic 
material described in the next section. The observations 
and results for both stars are summarized in Table [2j Be- 
low we provide a brief description of each of the instru- 
ments used, the data reduction, and the analysis proce- 
dure. 
We used the DS to conduct observations of both HAT- 



P-18 and HAT-P-19. This instrument delivers high- 
resolution spectra (A/AA ~ 35,000) over a single order 
centered on the Mg lb triplet (^5187 A). We measure the 
RV and stellar atmospheric param eters from the spect ra 
following the method described bv lTorres et al.l (J2002). 

We used FIES to conduct observations of HAT-P-19. 
We used the medium and the high-resolution fibers with 
resolving powers of A/AA « 46,000 and 67,000 respec- 
tively, giving a wavelength coverage of ~ 3600- 7400 A. 
The spectra were extracted and analyzed to measure the 
RV and stellar atmos pheric parameters follow ing the pro- 
cedures described by iBuchhave et ail (|2010[ ). The veloc- 
ities were corrected to the same system as the DS ob- 
servations (heliocentric velocities with the gravitational 
redshift of the Sun subtracted) using 10 observations of 
the velocity standard HD 182488 obtained on the same 
nights as observations of HAT-P-19. 

A single observation of HAT-P-19 was obtained with 
TRES. We used the medium-resolution fiber to obtain 
a spectrum with a resolution of A/AA w 44, 000 and a 
wavelength coverage of ~ 3900-8900 A. The spectrum was 
extracted and analyzed in a similar manner to the FIES 
observations. The velocity was corrected to the same 
system as the DS observations using a single TRES mea- 
surement of HD 182488 obtained on the same night. The 
velocity uncertainty reported in this case is our estimate 
of the systematic error based on the rms of multiple ob- 
servations of HD 182488 obtained on other nights close 
in time. 

Based on the reconnaissance spectroscopy observations 
we find that both systems have rms residuals consis- 
tent with no detectable RV variation within the preci- 
sion of the measurements. All spectra were single-lined, 
i.e., there is no evidence that either target star has a 
stellar companion. Additionally, both stars have sur- 
face gravity measurements which indicate that they are 
dwarfs. We note that for HAT-P-19 all three instruments 
yielded similar results for the RV and stellar parameters. 
There is a ~ 1 kms -1 difference between the DS and the 
TRES/FIES observations of HAT-P-19. The last DS ob- 
servation was obtained only four nights before the first 
FIES observation, and the DS and FIES data-sets each 
span significantly more than four nights, but do not show 
internal variations at the ~ lkms -1 level. We conclude 
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that the velocity difference between the instruments does 
not indicate a physical variation in the velocity of HAT- 
P-19. The DS observations are all weak, with only a few 
counts per pixel, and the sky velocity is slightly more 
negative than the system velocity in all cases. This ve- 
locity difference may be due to a systematic error in the 
DS velocities due to sky contamination. 

2.3. High resolution, high S/N spectroscopy 

We proceeded with the follow-up of each candidate by 
obtaining high-resolution, high-S/N spectra to character- 
ize the RV variations, and to refine the determination of 
the stellar parameters. These observations are summa- 
rized in Table [3J The RV measurements and uncertain- 
ties for HAT-P-18 are given in Table H and for HAT- 
P-19 in Table [5] The period- folded data, along with our 
best fit described below in Section |3l are displayed in 
Figure G2 for HAT-P-18, and in Figure H for HAT-P-19. 
For HAT-P-18, we exclude five RV measurements that 
are significant outliers from the best fit model. These 
points are all strongly affected by contamination from 
scattered moonlight (see Section 13.2.11) . Below we briefly 
describe the instruments used, the data reduction, and 
the analysis procedure. 

Observations were made of HAT-P-18 and H AT-P-19 
with the HIRES instrument ()Vogt et al.l 119941 ) on the 
Keck I telescope located on Mauna Kea, Hawaii. The 
width of the spectrometer slit was 0'.'86, resulting in a 
resolving power of A/AA ~ 55,000, with a wavelength 
coverage of ~3800-8000A. Exposures were obtained 
through an iodine gas absorption cell, which was used 
to superimpose a dense forest of I2 lines on the stellar 
spec trum and e stablish an accurate wavelength fiducial 
(see iMarcv fc B utler 1992). For each target two addi- 
tional exposures were taken without the iodine cell; in 
both cases we used the second, higher S/N, observation 
as the template in the reductions. Relative RVs in the 
sol ar system bar y centri c frame were derived as described 
by iButler et al.l ([1996), incorporating full modeling of 
the spatial and temporal variations of the instrumental 
profile. 

We also made use of the H igh-Dispersion Spectrograph 
(HDS; INoguchi et al.|[2002T) on the Subaru telescope on 
Mauna Kea, Hawaii to obtain high-S/N spectroscopic 
observations of HAT-P-19 from which we derived high- 
precision RV measurements. Observations were made 
over three consecutive nights using a slit width of 0'.'6, 
yielding a resolving power of A/AA ss 60,000. We 
used the I2b setup which provides a wavelength cover- 
age of ~ 3500 — 6200 A. To reduce the effect of changes 
in the barycentric velocity correction during an expo- 
sure, we limited exposure times to 15 minutes. As for 
Keck/HIRES, we made use of an iodine gas absorption 
cell to establish an accurate wavelength fiducial for each 
exposure. We also obtained six spectra without the io- 
dine cell, which were combined to form the template ob- 
servation. The spectra were extracted, and reduced to 
relative RVs in the solar syste m barycentr i c frame follo w- 
ing the methods described bv lSato et all (J20021 120051 ). 

In each figure we show also the relative S index, which 
is a measure of the chromospheric activity of the star 
derived from the flux in the cores of the Ca II H and K 
lines. This index was computed following the prescrip- 
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Fig. 3. — Top panel: Keck/HIRES RV measurements for 
HAT-P-18 shown as a function of orbital phase, along with our 
best-fit model (see Table [9j . Zero phase corresponds to the time 
of mid-transit. The center-of-mass velocity has been subtracted. 
Second panel: Velocity O—C residuals from the best fit. The error 
bars include a component from velocity jitt er (S .Oms -1 ) added in 
quadrature to the formal errors (see Section 13.311 . Third panel: Bi- 
sector spans (BS), with the mean value subtracted. The m easure - 
ment from the template spectrum is included (see Section l3.2.1|l . 
These measurements have not been corrected for contamination 
from moonlight; the corrected BS are shown in Figure [8] Bottom 
panel: Relative chromospheric activity index S measured from the 
Keck spectra. The formal errors on S based on photon statistics 
are comparable to the size of the displayed symbols. The scatter, 
however, is likely dominated by systematic errors in the measure- 
ments. Note the different vertical scales of the panels. 



tion given by iVaughan. Preston fc Wilson! (|1978P i , after 
matching each spectrum to a reference spectrum using 
a transformation that includes a wavelength shift and a 
flux scaling that is a polynomial as a function of wave- 
length. The transformation was determined on regions of 
the spectra that are not used in computing this indicator. 
Note that our relative S index has not been calib rated 
to the scale of IVaughan. Preston fc Wilsonl (|1978fl . We 
do not detect any significant variation of the index cor- 
related with orbital phase; such a correlation might have 
indicated that the RV variations could be due to stellar 
activity, casting doubt on the planetary nature of the 
candidate. 

2.4. Photometric follow-up observations 

In order to permit a more accurate modeling of the 
light curves, we conducted additional photometric ob- 
servations with the KeplerCam CCD camera on the 
FLWO 1.2 m telescope. The observations for each tar- 
get are summarized in Table [TJ 

The reduction of these images, including basic cali- 
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TABLE 2 
Summary of reconnaissance spectroscopy observations 



Instrument 



Date(s) 



Number of Spectra 



(K) 



log g* v sin i 

(cgs) (kms -1 ) 



7RV 
(kms -1 ) 



HAT-P- 


18 






DS 




2007 Sep- 


-2008 Mar 


HAT-P- 


19 






DS 

TRES 

FIES 




2008 Dec 

2009 Sep 
2009 Jan- 


-2009 Jan 

04 

-2009 Oct 



1 


4750 ± 100 


4.50 ±0.25 


1.5 ±1.0 


-11.92 ±0.28 


3 

1 
7 


4875 ± 125 
5000 ± 125 
4875 ± 133 


4.25 ±0.25 

4.5 ±0.25 

4.25 ±0.27 


5±5 

2±2 

2.1 ±2 


-21.2 ±0.5 
-20.20 ±0.05 
-20.22 ±0.02 



a The mean heliocentric RV of the target. 

TABLE 3 

Summary of high-resolution/high-SN spectroscopic 

observations 



Instrument 



Date(s) 



Number of 
RV obs. 



HAT-P-18 

Kcck/HIRES 
HAT-P-19 

Keck/HIRES 
Subaru/HDS 



2007 Oct-2010 Mar 

2009 Oct-2010 Feb 
2009 Aug 8-2009 Aug 10 



29 a 



13 

26 



a This number includes five outlier RV points which were ex- 
cluded from the analysis for HAT-P-18. 

bration, astrometry, and aperture photomet ry, was per- 
formed as described by iBakos et al.l (|2010( l We per- 
formed EPD and TFA to remove trends simultaneously 
with the light curve modeling; (for more details, see 
Section [3] and Bak os et"aLll2010f ). The final time series, 
together with our best-fit transit light curve model, are 
shown in the top portion of Figures [S] and [B] for HAT-P- 
18 and HAT-P-19 respectively; the individual measure- 
ments are reported in Tables [6] and [7] 

3. ANALYSIS 

3.1. Properties of the parent star 

Fundamental parameters for each of the host stars, 
including the mass (M*) and radius (i?*), which are 
needed to infer the planetary properties, depend strongly 
on other stellar quantities that can be derived spec- 
troscopically. For this we have relied on our template 
spectra obtained with the Keck/HIRES instrument, and 
the analysis packa ge known as S pectroscopy Made Easy 
(SME: IValenti fc Piskunov 1996), along with the atomic 
line database of IValenti fc Fischerl ((2005). For each star, 
SME yielded the following initial values and uncertain- 
ties (which we have conservatively increased by a factor 
of two to include our estimates of the systematic errors): 

• HAT-P-18 - effective temperature T off * =4850 ± 
75 K, stellar surface gravity log g± =4.7 ± 0.1 (cgs), 
metallicity [Fe/H] =+0.12±0.05dex, and projected 
rotational velocity vsmi =1.2 ± 0.5kms _1 . 

• HAT-P-19 - effective temperature T off * =5037 ± 
44K, stellar surface gravity logg* =4.7 ± 0.1 (cgs), 
metallicity [Fe/H] =+0.24±0.03dex, and projected 
rotational velocity vsmi =2.4 ± 0.5kms _1 . 



As discussed in Section 13. 2. 1| contamination from scat- 
tered moonlight affects the bisector spans and radial ve- 
locities measured for HAT-P-18 and the bisector spans 
measured for HAT-P-19. For HAT-P-18 the moon was 
below the horizon when the template used for the SME 
analysis was obtained, so it is not affected by contami- 
nation. For HAT-P-19 we estimate that scattered moon- 
light may have contributed ~ 0.1% of the total flux to 
the template spectrum used for SME analysis. The error 
in the parameters that results from this contamination is 
likely dwarfed by other systematic errors in the parame- 
ter determination. 

In principle the effective temperature and metallicity, 
along with the surface gravity taken as a luminosity in- 
dicator, could be used as constraints to infer the stellar 
mass and radius by comparison with stellar evolution 
models. However, the effect of log g* on the spectral line 
shapes is rather subtle, and as a result it is typically dif- 
ficult to determine accurately, so that it is a rather poor 
luminosity indicator in practice. For planetary transits a 
stronger constraint is often provided by the a/R+ normal- 
ized semimajor axis, which is closely related to p*, the 
mean stellar density. The quantity a/R* can be derived 
direc tly from the transi t lig ht curves (see iSozzetti et al.l 
120071 and also Section l3~3|) . This, in turn, allows us 
to improve on the determination of the spectroscopic 
parameters by supplying an indirect constraint on the 
weakly determined spectroscopic value of log g*, which 
removes degeneracies. We take this approach here, as de- 
scribed below. The validity of our assumption, namely 
that the best physical model describing our data is a 
planetary transit (as opposed to a blend), is shown later 
in Section 13.2.11 

For each system, our initial values of T er r*, log 5*, 
and [Fe/H] were used to determine auxiliary quantities 
needed in the global modeling of the follow-up photome- 
try and radial velocities (specifically, the limb-darkening 
coefficients). This modeling, the details of which are 
described in Section I3~3l uses a Monte Carlo approach 
to deliver the numerical probability distribution of a/R+ 
and other fit ted v ariable s. For further details we refer 
the reader to iPall (|2009bD . When combining a/R* (used 
as a proxy for luminosity) with assumed Gaussian dis- 
tributions for T c g+ and [Fe/H] based on the SME de- 
terminations, a comparison with stellar evolution mod- 
els allows the probability distributions of other stellar 
properties to be inferred, including logg*. Here we use 
the stellar evo lution calculations from Yonsei-Yale (YY; 
lYi et al.l 12001) for both stars. The comparison against 
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TABLE 4 

Relative radial velocities, bisector spans, and activity index 

measurements of hat-p-18. 



BJD a 


RV b 


°-RV C 


BS 


CBS 


S d 


TS 


(2,454,000+ 


(ms- 1 ) 


(ms -1 ) 


(ms -1 ) 


(ms -1 ) 






397.73420 . 


-11.23 


3.79 


-28.29 


13.90 


0.8595 


0.0099 


548.08129 . 






14.12 


20.54 


0.9267 


0.0051 


548.09633 . 


25.83 


2.26 


-2.20 


26.16 


0.9344 


0.0057 


549.07743 . 


22.80 


2.27 


3.36 


14.43 


0.9748 


0.0048 


602.83762 . 


9.95 


2.25 


67.02 


28.63 


1.0102 


0.0048 


602.99835 . 


15.96 


2.29 


12.79 


12.66 


1.0183 


0.0045 


603.83449 . 


20.24 


2.72 


141.82 


43.52 


1.0057 


0.0061 


604.08627 . 


10.30 


2.51 


-27.50 


9.88 


1.0266 


0.0049 


633.97708 . 


-21.65 


2.45 


9.96 


12.03 


1.0472 


0.0146 


635.98889 . 


20.26 


2.69 


52.95 


10.05 


0.9688 


0.0058 


639.02280 . 


-25.53 


2.29 


-8.96 


6.28 


0.9811 


0.0058 


641.97905 . 


31.77 


2.88 


-20.46 


13.72 


1.0007 


0.0075 


724.84068° 


6.49 


3.39 


-198.81 


36.35 


0.9206 


0.0158 


726.79620 . 


-27.63 


2.46 


-72.66 


14.74 


0.9544 


0.0091 


727.79088 . 


-26.34 


2.31 


-47.93 


16.52 


0.9805 


0.0078 


777.69605 . 


-16.44 


2.46 


26.52 


7.61 


0.9408 


0.0076 


778.69476 e 


35.06 


3.21 


45.44 


11.89 


0.8932 


0.0127 


779.70147" 


56.60 


4.98 


60.37 


20.83 


0.8978 


0.0212 


779.74051° 


61.15 


5.71 


134.75 


51.14 


0.8017 


0.0458 


865.15876 . 


-25.60 


4.88 


-31.22 


22.33 


1.0424 


0.0166 


955.02444 . 


14.38 


2.41 


-19.22 


14.02 


1.0448 


0.0056 


955.96262 . 


30.73 


2.64 


-17.35 


10.67 


1.0609 


0.0050 


964.09646 . 


-33.57 


2.32 


18.65 


11.88 


0.9697 


0.0047 


987.98665 . 


20.65 


2.88 


76.71 


20.29 


1.0195 


0.0057 


988.90490 . 


29.18 


2.82 


91.83 


24.70 


1.0630 


0.0058 


1109.76458° 


-49.58 


3.25 


-135.43 


28.38 


0.7973 


0.0308 


1016.99998 


30.57 


2.20 


-25.04 


8.82 


1.0197 


0.0076 


1041.96078 


-13.06 


3.12 


-12.39 


11.79 


1.0428 


0.0162 


1252.06235 


-3.90 


2.92 


-45.97 


12.38 


1.0616 


0.0065 


1252.07745 






-34.63 


6.47 


1.0437 


0.0055 


1261.11905 


-19.39 


2.39 


-28.25 


12.16 


1.0019 


0.0053 



Note. — Note that for the iodine-free template exposures we do not measure 
the RV but do measure the BS and S index. Such template exposures can be 
distinguished by the missing RV value. 

a Baryccntric Julian dates throughout the paper are calculated from Coordinated 
Universal Time (UTC) 

The zero-point of these velocities is arbitrary. An overall offset 7 rc i fitted to 
these velocities in Section l3.3l has not been subtracted. 

c Internal errors excluding the component of velocity jitter considered in 
Scction l3.3l 

Relative chromospheric activity index, not calibrated to the scale of 
Vaug han. Preston fc Wilson! U978I) . 
e Outlier RV measurements excluded from the analysis. 



the model isochrones was carried out for each of 10,000 
Monte Carlo trial sets for HAT-P-18, and 20,000 Monte 
Carlo trial sets for HAT-P-19 (see Section |3~3")) . Param- 
eter combinations corresponding to unphysical locations 
in the H-R diagram (41% of the trials for HAT-P-18 and 
31% of the trials for HAT-P-19) were ignored, and re- 
placed with another randomly drawn parameter set. For 
each system we carried out a second SME iteration in 
which we adopted the value of log g* so determined and 
held it fixed in a new SME analysis (coupled with a new 
global modeling of the RV and light curves), adjusting 
only T e ff^, [Fe/H], and wsini. This gave 

• HAT-P-18-T e B* =4803±80K, logg* =4.56±0.06, 
[Fe/H] =+0.10±0.08, andusini =0.5±0.5kms -1 . 

• HAT-P-19-T cS + =4990±130K, log 5 * =4.53±0.06 
(fixed), [Fe/H] =+0.23 ± 0.08, and vsini =0.7 ± 
0.5 kms -1 . 

In each case the conservative uncertainties for T e ff* and 
[Fe/H] have been increased by a factor of two over their 



formal values, as before. For each system, a further iter- 
ation did not change log g± significantly, so we adopted 
the values stated above, together with the new logg* 
values resulting from the global modeling, as the final 
atmospheric properties of the stars. They are collected 
in Table [8] for both stars. 

With the adopted spectroscopic parameters the model 
isochrones yield the stellar mass and radius, and other 
properties. These are listed for each of the systems in Ta- 
ble [SJ According to these models HAT-P-18 is a dwarf 
star with an estimated age of 12.4t 64 Gyr, and HAT-P- 
19 is a dwarf star with an estimated age of 8.8 ± 5.2 Gyr. 
The inferred location of each star in a diagram of a/ R* 
versus T e gf*, analogous to the classical H-R diagram, is 
shown in Figure [7] The stellar properties and their la 
and 2<t confidence ellipsoids are displayed against the 
backdrop of model isochrones for a range of ages, and the 
appropriate stellar metallicity. For comparison, the loca- 
tions implied by the initial SME results are also shown 
with triangles. 
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TABLE 5 

Relative radial velocities, bisector spans, and activity index 

measurements of hat-p-19. 



BJD a 


RV b 


<mv c 


BS 


CBS 


S d 


cs 


Inst. 


(2,454,000+ 


(ms -1 ) 


(ms- 1 ) 


(ms -1 ) 


(ms" 1 ) 








1052.04317 


-83.78 


7.01 


7.91 


33.17 


1.042 


0.324 


Subaru 


1052.05438 


-81.45 


6.85 


9.71 


30.12 


0.905 


0.294 


Subaru 


1052.06559 


-71.03 


6.92 


14.66 


27.93 


1.016 


0.315 


Subaru 


1052.07680 


-92.65 


6.96 


-14.95 


29.17 


1.109 


0.328 


Subaru 


1052.08801 


-76.43 


6.91 


-8.57 


29.23 


1.059 


0.330 


Subaru 


1052.09923 


-72.59 


6.90 


-18.67 


31.68 


1.079 


0.329 


Subaru 


1052.11044 


-83.27 


6.87 


-11.72 


33.40 


1.006 


0.304 


Subaru 


1052.12165 


-87.81 


6.88 


-8.79 


35.31 


1.108 


0.340 


Subaru 


1052.99045 






-25.11 


34.01 


0.982 


0.313 


Subaru 


1053.00513 






-15.75 


37.37 


1.018 


0.318 


Subaru 


1053.01981 






-38.58 


39.13 


1.083 


0.330 


Subaru 


1053.03449 






-6.01 


32.64 


1.028 


0.332 


Subaru 


1053.04918 






19.45 


23.05 


0.986 


0.315 


Subaru 


1053.06386 






-22.72 


35.09 


1.126 


0.347 


Subaru 


1053.11425 


-65.38 


6.90 


15.97 


38.62 


0.962 


0.309 


Subaru 


1053.12894 


-79.28 


6.85 


13.01 


37.16 


0.948 


0.313 


Subaru 


1053.14023 


-67.57 


7.11 


37.01 


30.50 


1.011 


0.331 


Subaru 


1053.93729 


-9.14 


6.91 


-11.41 


35.52 


1.111 


0.338 


Subaru 


1053.94850 


0.07 


6.93 


2.85 


33.24 


1.023 


0.320 


Subaru 


1053.95970 


-21.48 


6.88 


3.11 


36.49 


1.085 


0.334 


Subaru 


1053.97091 


-14.40 


6.87 


-1.28 


34.68 


1.040 


0.314 


Subaru 


1053.98213 


-32.56 


6.87 


12.91 


33.85 


0.968 


0.301 


Subaru 


1053.99335 


-10.54 


6.84 


-12.62 


38.43 


0.926 


0.304 


Subaru 


1054.00457 


-24.54 


6.83 


-10.77 


32.87 


0.933 


0.306 


Subaru 


1054.04057 


-23.88 


6.80 


5.45 


35.97 


0.855 


0.273 


Subaru 


1054.05178 


-14.33 


6.78 


35.76 


34.71 


0.907 


0.295 


Subaru 


1054.06300 


-13.93 


6.82 


32.75 


23.55 


0.870 


0.285 


Subaru 


1054.07421 


-20.91 


6.81 


-4.30 


35.97 


1.009 


0.323 


Subaru 


1054.08543 


-9.03 


6.81 


-8.63 


35.77 


0.947 


0.310 


Subaru 


1054.09665 


-18.96 


6.75 


-3.25 


35.11 


0.953 


0.304 


Subaru 


1054.10787 


-17.49 


6.78 


5.49 


35.24 


0.967 


0.308 


Subaru 


1054.11908 


-14.29 


6.82 


7.07 


34.89 


0.937 


0.297 


Subaru 


1107.06075 






7.74 


7.81 


1.048 


0.011 


Keck 


1107.07559 


4.11 


2.77 


18.29 


9.69 


1.053 


0.044 


Keck 


1108.99050 


-69.19 


3.09 


54.83 


18.48 


1.009 


0.014 


Keck 


1112.11372 


-51.05 


3.13 


35.61 


17.65 


0.928 


0.020 


Keck 


1134.04526 


7.06 


2.60 


-33.21 


6.43 


1.013 


0.024 


Keck 


1136.01174 


-31.79 


3.10 


169.45 


34.81 


0.928 


0.034 


Keck 


1172.81382 


-61.98 


2.36 


-14.18 


3.85 


1.075 


0.011 


Keck 


1188.84358 


-27.51 


2.19 


-22.96 


5.20 


1.014 


0.010 


Keck 


1190.80886 


53.47 


2.62 


-32.46 


11.07 


0.998 


0.011 


Keck 


1192.90496 


-36.33 


2.34 


-92.65 


14.20 


0.660 


0.020 


Keck 


1193.80419 


28.07 


2.18 


-20.27 


8.65 


1.101 


0.011 


Keck 


1198.81373 


59.29 


2.49 


-50.27 


11.66 


1.020 


0.010 


Keck 


1250.72949 






-14.06 


11.36 


1.072 


0.011 


Keck 


1250.74616 


70.58 


2.56 


-4.56 


10.41 


1.025 


0.015 


Keck 


1251.72400 


44.90 


2.60 


-1.30 


10.91 


1.056 


0.016 


Keck 



Note. — Note that for the iodine-free template exposures we do not measure the RV 
but do measure the BS and S index. Such template exposures can be distinguished by 
the missing RV value. 

a Baryccntric Julian dates throughout the paper are calculated from Coordinated Uni- 
versal Time (UTC) 

The zero-point of these velocities is arbitrary. An overall offset 7 rc i fitted to these 
velocities in Scction l3.3l has not been subtracted. 
c Internal errors excluding the component of velocity jitter considered in Scction l3.3l 

Relative chromospheric activity index, not calibrated to the scale of 
Vaughan, Preston & Wilson (197^). Note the values for the Keck and Subaru obser- 
vations have independently been scaled to have a mean of 1.0. 



The stellar evolution modeling provides color indices 
that may be compared against the measured values as a 
sanity check. For each star, the best available measure- 
ments are the near-infrared magnitudes from the 2MASS 
Catalogue (jSkr utskic ct al. 2006), which are given in 
Table [3] These are converted to the photometric sys- 
tem of t he models ( ESQ system) using the transforma- 
tions by iCarpenterl (|2001l ). The resulting color index is 
J -K = 0.623 ± 0.036, and J - K = 0.583 ± 0.031 for 
HAT-P-18, and HAT-P-19 respectively. These are both 



within ler of the predicted values from the isochrones of 
J - K = 0.61 ± 0.02 and J - K = 0.57 ± 0.03. The dis- 
tance to each object may be computed from the absolute 
K magnitude from the models and the 2MASS K s mag- 
nitudes, which has the advantage of being less affected 
by extinction than optical magnitudes. The results are 
given in Table [51 where in each case the uncertainty ex- 
cludes possible systematics in the model isochrones that 
are difficult to quantify. 
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Fig. 4. — Top panel: High-precision RV measurements for 
HAT-P-19 from Keck/HIRES (dark filled circles) and Subaru/HDS 
(light filled squares) shown as a function of orbital phase, along 
with our best- fit model (see Table [5). Zero phase corresponds to 
the time of mid- transit. The center-of-mass velocity and a linear 
trend (second panel) have been subtracted. Second panel: Velocity 
O — C residuals from the best fit single Keplerian orbit model as a 
function of time. The residuals show a linear trend, indicative of a 
third body in the system. Note that the velocity zero points of the 
Subaru and Keck observations are independently free parameters. 
The linear trend is thus not constrained by the Subaru observations 
which span only 3 days. Third panel: Velocity O—C residuals from 
the best fit including both the Keplerian orbit and linear trend, 
shown as a function of orbital phase. The error bars include a 
component from velocity jitter (6.7 ms" 1 ) added in quadrature to 
the formal errors (see Section 13.311 . Fourth panel: Bisector spans 
(BS), with the mean value subtracted. The meas urement from the 
template spectrum is included (see Section I3.2.1J I. These measure- 
ments have not been corrected for contamination from moonlight; 
the corrected BS are shown in Figure [8] Bottom panel: Relative 
chromospheric activity index S measured from the Keck spectra. 
The formal errors on S based on photon statistics are comparable 
to the size of the displayed symbols. The scatter, however, is likely 
dominated by systematic errors in the measurements. Note the 
different vertical scales of the panels. Observations shown twice 
are represented with open symbols. 
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Fig. 5. — Unbinned transit light curves for HAT-P-18, acquired 
with KeplerCam at the FLWO 1.2 m telescope. The li ght c urves 
have been EPD and TFA processed, as described in § 13.31 The 
dates of the events are indicated. The second curve is displaced 
vertically for clari ty. Our best fit from the global modeling de- 
scribed in Section 13.31 is shown by the solid lines. Residuals from 
the fits are displayed at the bottom, in the same order as the top 
curves. The error bars represent the photon and background shot 
noise, plus the readout noise. 
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Fig. 6. — Similar to Figurc[5] here we show the follow-up light 
curves for HAT-P-19. 



3.2. Rejecting Blend Scenarios 

Our initial spectroscopic analyses discussed in 
Section [2~2l and Section [231 rule out the most obvious as- 
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Fig. 7.— Left: model isochrones fromlYT ct al. (2001]) for the measured metallicity of HAT-P-18, [Fe/H]= +0.10, and ages of 0.2, 0.5, and 
1 to 14Gyr, in 1 Gyr increments (left to right). The adopted values of T e g* and a/Ri, are shown together with their 1<t and 2<r confidence 
ellipsoids. The initial values of T e ff+ and a/Ft* from the first SME and light curve analyses are represented with a triangle. Right: same 
as left, here we show the results for HAT-P-19, with [Fe/H]= +0.23, and ages of 1 to 13 Gyr in steps of 1 Gyr (left to right). 



TABLE 6 
High-precision differential photometry of HAT-P-18. 



BJD 


Mag a 


""Mag 


Mag(orig) b 


Filter 


(2,400,000+) 










54582.72249 


-0.00198 


0.00093 


10.79530 


i 




54582.72334 


0.00006 


0.00093 


10.79800 






54582.72420 


-0.00247 


0.00093 


10.79660 


i 




54582.72506 


0.00024 


0.00093 


10.79950 


i 




54582.72592 


0.00344 


0.00093 


10.80160 






54582.72677 


-0.00104 


0.00093 


10.79740 


i 




54582.72760 


0.00148 


0.00092 


10.79980 


i 




54582.72848 


0.00027 


0.00092 


10.79820 






54582.72933 


0.00002 


0.00092 


10.79590 


i 




54582.73102 


-0.00066 


0.00092 


10.79560 


i 



Note. — This tabic is available in a machine-readable form in the 
online journal. A portion is shown here for guidance regarding its form 
and content. 

a The out-of-transit level has been subtracted. These magnitudes have 
been subjected to the EPD and TFA procedures, carried out simulta- 
neously with the transit fit. 

Raw magnitude values without application of the EPD and TFA 
procedures. 



TABLE 7 
High-precision differential photometry of HAT-P-19. 



BJD 


Mag a 


°"Mag 


Mag(orig) b 


Filter 


(2,400,000+) 










55111.59368 


0.02677 


0.00139 


11.66270 


? 




55111.59465 


0.02295 


0.00130 


11.65390 


i 




55111.59562 


0.02551 


0.00142 


11.65450 


i 




55111.59853 


0.02643 


0.00133 


11.66000 


1 




55111.59949 


0.02348 


0.00120 


11.65220 


i 




55111.60045 


0.02275 


0.00140 


11.65930 


1 




55111.60140 


0.02615 


0.00126 


11.65320 


1 




55111.60236 


0.02510 


0.00126 


11.65380 


i 




55111.60333 


0.02741 


0.00130 


11.66230 


i 




55111.60430 


0.02038 


0.00121 


11.64810 


i 



Note. — This table is available in a machine-readable form in 
the online journal. A portion is shown here for guidance regarding 
its form and content. 

a The out-of-transit level has been subtracted. These magnitudes 
have been subjected to the EPD and TFA procedures, carried out 
simultaneously with the transit fit. 

Raw magnitude values without application of the EPD and TFA 
procedures. 



trophysical false positive scenarios. However, more sub- 
tle phenomena such as blends (contamination by an un- 
resolved eclipsing binary, whether in the background or 
associated with the target) can still mimic both the pho- 
tometric and spectroscopic signatures we see. In the fol- 
lowing sections we consider and rule out the possibility 
that such scenarios may have caused the observed pho- 
tometric and spectroscopic features. 

3.2.1. Spectral line-bisector analysis 

Following iTorres et al.1 (|2007| ). we explored the possi- 
bility that the measured radial velocities are not real, 
but are instead caused by distortions in the spectral line 
profiles due to contamination from a nearby unresolved 
eclipsing binary. A bisector analysis for each system 
based on the Keck spectra (and the Subar u spectra for 
HAT- P-19) was done as described in §5 of iBakos et al.l 
(|200l . 

Each system shows excess scatter in the bisector spans 
(BS), above what is expected from the measurement er- 
rors (see Figure |3l third panel, and Figure IU fourth 
panel). For HAT-P-18 there may be a slight correla- 
tion between the RV and the BS, while for HAT-P-19 no 
correlation is apparent. Such a correlation could indi- 
cate that the photometric and spectroscopic signatures 
are due to a blend scenario rather than a single planet 
transiting a single star. We note that for HAT-P-19, the 
Keck spectra show excess BS variation, while the Subaru 
spectra do not. 

Following our ea rlier work ()Kovacs et al.l 120101 : 
lHartman et al.l 120091 ) we investigated the effect of con- 
tamination from moonlig h t on the measured BS val- 
ues. As in IKovacs et al.l (|2010t) . we estimate the ex- 
pected BS value for each spectrum by modeling the spec- 
trum cross-correlation function (CCF) as the sum of two 
Lorentzian functions, shifted by the known velocity dif- 
ference between the star and the moon, and scaled by 
th eir expected flux ratio (estimated following equation 3 
of lHartman et al.l 12009) . We refer to the simulated BS 
value as the sky contamination factor (SCF). We find 
a strong correlation between the SCF and BS for both 
systems (see Figure [5]). After correcting for this correla- 
tion, we find that the BS show no significant variations, 
and the correlation between the RV and BS variations 
is insignificant for both systems. Therefore, we conclude 
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TABLE 8 
Stellar parameters for HAT-P-18 and HAT-P-19 



Parameter 


HAT-P-18 


HAT-P-19 


Source 


Spectroscopic properties 






Trf* (K) 


4803 ± 80 


4990 ± 130 


SME a 


[Fe/H] 


+0.10 ±0.08 


±0.23 ±0.08 


SME 


iisini (kms 1 ) . . . 


0.5 ±0.5 


0.7 ±0.5 


SME 


%ic (kms -1 ) 


2.53 


2.81 


SME 


v mic (kms -1 ) 


0.85 


0.85 


SME 


7RV (kms -1 ) 


-11.92 ±0.28 


-20.22 ± 0.02 


DS/FIES b 


Photometric properties 








V (mag) 


12.759 


12.901 


TASS 


V-I c (mag) 


1.18 ±0.13 


1.04 ±0.13 


TASS 


J (mag) 


10.822 ±0.020 


11.095 ± 0.020 


2MASS 


ff (mag) 


10.340 ±0.019 


10.644 ± 0.022 


2MASS 


K s (mag) 


10.234 ±0.017 


10.546 ±0.019 


2MASS 


Derived properties 








M, (M ) 


0.770 ± 0.031 


0.842 ± 0.042 


YY+a/fc+SME c 


-R* (Rq) 


0.749 ±0.037 


0.820 ± 0.048 


YY+a/JL+SME 


log 9* (cgs) 


4.57 ±0.04 


4.54 ±0.05 


YY+a/R^+SME 


i* (i©) 


0.27 ±0.04 


n o 7 +0.08 
v -°' -0.06 


YY+a/fl*+SME 


My (mag) 


6.50 ±0.19 


6.08 ±0.24 


YY+a/i?*+SME 


M K (mag,ESO) . . 


4.17 ±0.12 


3.92 ±0.15 


YY+a/i?*+SME 


Age (Gyr) 


12.4±tt 


8.8 ±5.2 


YY+a/R*+SME 


Distance (pc) 


166 ±9 


215 ±15 


YY+a/R*+SME 



a SME — "Spectroscopy Made Easy" package for the analysis of high-resolution 
spectra (Valcnti & Piskunov 1996). These parameters rely primarily on SME, but 
have a small dependence also on the iterative analysis incorporating the isochronc 
search and global modeling of the data, as described in the text. 
b Based on DS observations for HAT-P-18 and FIE S observations for HAT-P-19. 
c YY+a/fl,+SME = Based on the YY isochrones jYi et al.l feOOl), a/R* as a lumi- 
nosity indicator, and the SME results. 



that the the velocity variations are real for both stars, 
and that both stars are orbited by close-in giant plan- 
ets. An independent method for arriving at this same 
conclusion is also presented in the following section. 

We have also investigated the effect of sky contami- 
nation on the measured RVs. The expected RV due to 
sky contamination for a given spectrum is estimated by 
finding the peak of the simulated CCF. Note that the 
real RV measurements are obtained by directly model- 
ing the spectra and not by performing cross-correlation. 
We therefore only expect a crude agreement between the 
estimated RVs due to sky contamination, and the real 
RVs. Figure |9] compares the expected RVs to the mea- 
sured RV residuals from the best-fit model for HAT-P-18 
and HAT-P-19. For HAT-P-18 we find a rough correla- 
tion between the estimated and measured RV residuals. 
Five of the RV measurements which are significant out- 
liers from the best-fit model are rejected. These spectra 
are also among the most strongly affected by sky con- 
tamination. For HAT-P-19 the values do not appear to 
be correlated. 

3.2.2. Blend Modeling of the Photometry 

As an independent test on the possibility that the 
observations for either HAT-P-18 or HA T-P-19 could 
be ca us ed by a blend s cenar io, we follow T orres et al.l 
(|2005t l. lHartman etHI (|2009f ). and iBakos et all (|2010f l 
in attempting to model the photometric observations for 
each object as either a hierarchical triple system, or a 
blend between a bright foreground star and a background 
eclipsing binary. We will show that for both HAT-P- 
18 and HAT-P-19 blend scenarios that do not include a 



transiting planet may be rejected from the photometric 
observations alone. We consider 5 possibilities: 



1. One star orbited by a planet, 

2. Hierarchical system, 3 stars, 
eclipsing, 



2 fainter stars are 



3. Hierarchical system, 2 stars, 1 planet, planet orbits 
the fainter star, 

4. Hierarchical system, 2 stars, 1 planet, planet orbits 
the brighter star, 

5. Chance alignment, 3 stars, 2 background stars are 
eclipsing. 

Here case 1 is the fiducial model to which we compare the 
various blend models. We model the observed follow-up 
and HATNet light curves (including only points that are 
within one transit duration of the primary transit or sec- 
ondary eclipse assuming zero eccentricity) together with 
the 2MASS and TASS photometry. In all cases we vary 
the distance to the brightest star in the system, parame- 
ters allowing for dilution in the HATNet light curves, and 
we include simultaneous EPD and TFA in fitting the light 
curves (see Section l3~3j) . We draw the ste llar radii and 
magn itudes from the Padova isochrones (Girar di et al. 
2000), extended below 0.15 M Q with the iBaraffe et aL 
(1998) isochrones. We use these rather than the YY 
isochrones for this analysis because of the need to allow 
for stars with M < 0.4 M©, which is the lower limit avail- 
able for the YY models. We use the JKTEBOP program 
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Fig. 8. — A: BS vs SCF for Keck spectra of HAT-P-18. The strong correlation between these quantities indicates that much of the BS 
variation for this object can be accounted for by changes in the sky contamination of the spectra. B: BS vs orbital phase for HAT-P-18 
after fitting and subtracting a linear relation between BS and SCF. The residual BS is uncorrelated with the orbital phase, justifying our 
conclusion that this system is not a blend. Bottom: BS vs SCF for Keck spectra of HAT-P-19. Again the BS variation for this target can 
be accounted for by the changing sky contamination. The Subaru spectra of HAT-P-19 do not show significant BS variations; these spectra 
were taken within a span of 3 days, under similar sky conditions. D: Same as panel B, here shown for HAT-P-19. 



(|Southworth et al1l2004al lbT) which is based on the Eclips- 
ing Binary Orbit Program (EBOP ; iPopper & Etze]|ll98H 
lEtzellll981b iNelson fc Davislll972f ) to generate the model 
light curves. We optimize the free parameters using the 
Downhill Simplex Algorithm together with the classical 
linear least squares algorithm for the EPD and TEA pa- 
rameters. We rescale the errors for each light curve such 
that x 2 P er degree of freedom is 1.0 for the out of tran- 
sit portion of the light curve. Note that this is done 
prior to applying EPD/TEA corrections for systematic 
errors. As a result, the \ 2 P er degree of freedom is less 
than 1.0 for many of the best-fit models discussed be- 
low. If the rescaling is not performed, the difference in 
X 2 between the best-fit models is even more significant 
than what is given below, and the blend models may 
be rejected with even higher confidence. For HAT-P- 
18 we fix the mass, age, and [Fe/H] metallicity of the 
brightest star in the system to 0.76 A/©, 12.4 Gyr, and 
+0.10 respectively to reproduce the effective tempera- 
ture, metallicity, and surface gravity of the bright star 
as determined from the SME analysis when using the 
Padova isochrones. For HAT-P-19 we fix the mass, age, 
and metallicity to 0.83 M@, 8.8 Gyr, and +0.23 respec- 
tively. 

Case 1: 1 star, 1 planet: In addition to the param- 
eters mentioned above, in this case we vary the radius 
of the planet and the impact parameter of the transit. 



For HAT-P-18 the best-fit model has x 



18,Cascl 



1542.4 



for 1712 degrees of freedom. For HAT-P-19 the best-fit 
model has xfg easel = 2710.0 for 2746 degrees of free- 
dom. The parameters that we obtain for both objects 
are comparable to those obtained from the global mod- 
elling described in Section 13.31 

Case 2: Hierarchical system, 3 stars: For case 2 we 
vary the masses of the eclipsing components, and the im- 
pact parameter of the eclipse. We take the radii and mag- 
nitudes of all three stars from the same isochrone. For 
HAT-P-18 we find Xi 8 ,Ca SC 2 = 1566 - 9 for ml degrees of 
freedom, while for HAT-P-19 we find X?g,Ca S e2 = 2753.6 
for 2745 degrees of freedom. The best-fit model for HAT- 
P-18 consists of a 0.76 M© star that is blended with a 
eclipsing binary with components of mass 0.74 M© and 
0.12 M Q . For HAT-P-19 the best-fit model consists of 
two equal 0.83 M© stars with a 0.13 M© M dwarf eclips- 
ing one of the two K stars. For both HAT-P-18 and HAT- 
P-19 the best-fit case 2 model has higher x 2 with fewer 
degrees of freedom than the best-fit case 1 model, so for 
both objects the case 1 model is prefered. To establish 
the significance at which we may reject the Case 2 model 
in favor of the Case 1 model, we conduct Monte Carlo 
simulations which account for the possibility of uncor- 
rcctc d systematic e rrors in the light curves as described 
in lHartman et al.l (|2009T >. For HAT-P-18 we reject the 
best-fit Case 2 model at the ~ 4.4er confidence level, 
while for HAT-P-19 we reject the best-fit Case 2 model 
at the ~ 5.1<7 confidence level. We also note that for 
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Fig. 9. — Top: Residual RV from best-fit model vs estimated 
RV due to sky contamination for Keck spectra of HAT-P-18. There 
is a rough correlation between the values. Open circles show RV 
outliers which were excluded when modeling the orbit, filled circles 
show all other measurements. The excluded measurements come 
from observations which are expected to be strongly contaminated 
by the sky. Bottom: Same as above, here we show the results 
for HAT-P-19. In this case there is no apparent correlation be- 
tween the observed RV residuals and the expected RVs due to sky 
contamination. 

both objects the only hierarchical triple stellar system 
that could potentially fit the photometric observations 
is a system where the two brightest stars have nearly 
equal masses. Because both HAT-P-18 and HAT-P-19 
have narrow spectral lines (i>sini =0.5 ± 0.5 km s^ 1 and 
0.7 ±0.5 km s -1 respectively), a second component with 
a luminosity ratio close to one and a RV semi-amplitude 
of several tens of km s _1 would have easily been detected 
in the spectra of these objects. 

Case 3: Hierarchical system, 2 stars, 1 planet, planet 
orbits fainter star. In this scenario the system contains 
a transiting planet, but it would have a radius that is 
larger than what we infer assuming there is only one star 
in the system. For this case we vary the mass of the faint 
planet-hosting star, the radius of the planet, and the im- 
pact parameter of the transit. We assume the mass of 
the planet is negligible relative to the mass of its faint 
host star. For HAT-P-18 the best-fit Case 3 model has 
Xi8,casc3 = 1590.0 while for HAT-P-19 the best-fit Case 3 
model has x?g Case3 = 2726.1. For both objects the best- 
fit is when the two stars in the system are of equal mass. 
Repeating the Monte Carlo simulations to determine the 
statistical significance of this x 2 difference, we find that 
for HAT-P-18 we may reject the Case 3 model in favor of 
the Case 1 model at the ~ 4.9<r confidence level, while for 
HAT-P-19 we may reject the Case 3 model at the ~ 3.6<r 



confidence level. As for the Case 2 model, the only Case 
3 models that could potentially fit the photometric ob- 
servations for either HAT-P-18 or HAT-P-19 are models 
where both stars in the system have nearly equal mass. 
The narrow spectral lines for both HAT-P-18 and HAT- 
P-19 means that the systemic velocities for the putative 
binary star companions would need to be very similar to 
those of the brighter stars not hosting the planets (within 
~ lkms -1 ) for the secondary stars to have gone unde- 
tected in any of our spectroscopic observations. 

Case 4- Hierarchical system, 2 stars, 1 planet, planet 
orbits brighter star. As in Case 3, in this scenario the 
system constains a transiting planet, but it would have a 
radius that is larger than what we infer assuming there 
is only one star in the system. For this case we vary the 
mass of the faint contaminating star, the radius of the 
planet, and the impact parameter of the transit. Again 
we assume the mass of the planet is negligible relative 
to the mass of its host star. For HAT-P-18 the best-fit 
model occurs when the mass of the contaminating star is 
negligible with respect to the mass of the planet-hosting 
star (which effectively corresponds to the Case 1 model) , 
with x?8 Casc4 increasing as the mass of the faint com- 
panion is increased. We find that a fainter companion 
with M > 0.62 M & is rejected at the 3a confidence level, 
while a companion with M > 0.51 M© is rejected at the 
2cr confidence level. This corresponds to 3cr and 2tr up- 
per limits on the V-band luminosity ratios of a possible 
contaminating star of 0.25 and 0.09 respectively. At the 
2cr level, the radius of HAT-P-18b could be 5% larger 
than what we measure in Section 13.31 if there is an unde- 
tected faint companion star. For HAT-P-19 we find that 
a fainter companion with M > 0.81 M© is rejected at the 
3(7 level, while a fainter companion with M > 0.69 M© is 
rejected at the 2cr level. This corresponds to 3a and 2a 
upper limits on the V-band luminosity ratios of a possi- 
ble contaminating star of 0.84 and 0.27 respectively. At 
the 2cr level, the radius of HAT-P-19b could be ~ 18% 
larger than what we measure in Section 13.31 if there is 
an undetected faint companion star. While \ 2 generally 
increases when the mass of the faint star increases, its 
minimum actually occurs when the faint star has a mass 
of ~ 0.4 — 0.5 M©, where the value of x 2 is 4 less than the 
value when the faint companion is excluded (the Case 1 
model). The difference is too low to be statistically sig- 
nificant, but it is nonetheless interesting that this object 
also exhibits a linear drift in its radial velocity, which 
might indicate the presence of a low-mass stellar com- 
panion. If there is a 0.45 M© faint companion, the ra- 
dius of HAT-P-19b would be ~ 3% larger than what we 
measure in Section 13.31 

Case 5: Chance alignment, 3 stars, background stars 
are eclipsing: For case 5 we vary the masses of the two 
eclipsing stars, the impact parameter of their eclipses, 
the age of the background system, the metallicity of the 
background system, and the difference in distance mod- 
ulus between the foreground star and the background bi- 
nary AV. For HAT-P-18 we find that the best-fit model 
has Xi8 Caso5 = 1561.2 and consists of a "background" 
binary at AV = with a primary component that has 
the same mass as the foreground star. This is effectively 
Case 2, except that the age and metallicity of the bi- 
nary are allowed to vary, so that the result has a slightly 
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lower x 2 value than the best-fit Case 2 model. The value 
of xf 8 Q ase5 steadily increases with AV. For HAT-P-18 
the best-fit Case 5 model may be rejected in favor of the 
Case 1 fiducial model at the ~ 4.1<r confidence level. We 
may reject models with AV > 0.2 with greater than 5a 
confidence. The 5a lower limit on the V-band luminosity 
ratio between the primary component of the background 
binary and the foreground star is ~ 0.75. Such a system 
would have easily been identified and rejected as a spec- 
troscopic double-lined object in either the Keck or TRES 
spectra. For HAT-P-19 we find that the best-fit model 
has X19 Cases — 2727.7 and consists of a background bi- 
nary at AV = 0.1 mag. The primary component of the 
binary has a mass of 0.85 M Q , while the secondary has a 
mass of 0. 14 M Q . The binary system has an age of 9 Gyr 
and a metallicity of [Fe/H]= +0.4. This model can be 
rejected in favor of the fiducial model at the 3.5a con- 
fidence level. We note that we can reject models with 
AV > 0.9 mag at the 5a level. The 5a lower limit on 
the V-band luminosity ratio between the primary compo- 
nent of the background binary and the foreground star is 
~ 0.58. As for HAT-P-18, such a system would have eas- 
ily been identified and rejected as a spectroscopic double- 
lined object in either the Keck or TRES spectra. We con- 
clude that for HAT-P-18 and HAT-P-19 a blend model 
consisting of a single star and a background eclipsing bi- 
nary is inconsistent with the photometric observations at 
the 3-4cr level, and any models that are not inconsistent 
with greater than 5a confidence are inconsistent with the 
spectroscopic observations. This reaffirms our conclusion 
in the previous section about the true planetary nature 
of the signals in both HAT-P-18 and HAT-P-19. 

3.3. Global modeling of the data 

This section describes the procedure we followed for 
each system to model the HATNet photometry, the 
follow-up photometry, and the radial velocities simulta- 
neously. Our model for the follow-up light curves used 
analytic formulae based on iMandel fc Agoll (|2002fi for 
the eclipse of a star by a planet, with limb darkening 
being prescribed by a quadratic law. The limb darken- 
ing coefficients for the Sloan i-band and Sloan q-b and 
were interpolated from the tables by IClaretl ()2004f ) for 
the spectroscopic parameters of each star as determined 
from the SME analysis ( Section 13. II) . The transit shape 
was parametrized by the normalized planetary radius 
p = Rp/B-t, the square of the impact parameter b 2 , and 
the reciprocal of the half duration of the transit C/-R*- 
We chose these parameters because of their simple ge- 
ometric meanings a nd the fact that t hese show negligi- 
ble correlations (see lBakos et al.ll2010D . The relation be- 
tween C/-R* and the quantity a/R±, used in Section [37X1 
is given by 



a/R+ = P/27r(C/-R*)Vl - &Vl - e 2 /(l + esinw) (1) 

(see, e.g., iTinglev fc Sackettl 120051) . Our model for 
the HATNet data w as a simplified version of the 
IMandel fe Agoll (|2002t ) analytic functions (an expansion 
in terms o f Legendre p o lynom ials), for the reasons de- 
scribed in IBakos et al.l (I2010D . Following the formal- 
ism presented by [Pal (2009), the RVs were fitted with 
an eccentric Keplerian model parametrized by the semi- 
amplitude K and Lagrangian elements k = e cos u and 



h = esinw, in which ui is the longitude of periastron. 

We assumed that there is a strict periodicity in the in- 
dividual transit times. For each system we assigned the 
transit number Ntr — to the first complete follow-up 
light curve. For HAT-P-18b this was the light curve 
obtained on 2008 Apr 25, and for HAT-P-19b this was 
the light curve obtained on 2009 Dec 1. The adjustable 
parameters in the fit that determine the ephemeris were 
chosen to be the time of the first transit center observed 
with HATNet (T c _ 7 i, and T c _ 2 04 for HAT-P-18b, and 
HAT-P-19b respectively) and that of the last transit cen- 
ter observed with the FLWO 1.2 m telescope (T c + 69, and 
T c , for HAT-P-18b, and HAT-P-19b respectively). We 
used these as opposed to period and reference epoch in 
order to minim ize correlations between parameters (see 
iPal et aD 120081 ). Times of mid-transit for intermediate 
events were interpolated using these two epochs and the 
corresponding transit number of each event, N tr - For 
HAT-P-18b, the eight main parameters describing the 
physical model were thus the times of first and last 
transit center, R p /R+, b 2 , (,/R*, K, k = ecosw, and 
h = esinw. For HAT-P-19b we included as a ninth pa- 
rameter a velocity acceleration term to account for an 
apparent linear drift in the velocity residuals after fit- 
ting for a Keplerian orbit. Three additional parameters 
were included for HAT-P-18b that have to do with the in- 
strumental configuration. For HAT-P-19b, six additional 
parameters were included. These include the HATNet 
blend factors B- mst (one for each HATNet field for HAT- 
P-19b), which accounts for possible dilution of the transit 
in the HATNet light curve from background stars due to 
the broad PSF (20" FWHM), the HATNet out-of-transit 
magnitude M^HATNct (also one for each HATNet field for 
HAT-P-19b), and the relative zero-point 7 rc i of the Keck 
RVs (and the Subaru RVs for HAT-P-19b). 

We extended our physical model with an instrumental 
model that describes brightness variations caused by sys- 
tematic errors in the measurements. This was done in a 
similar fashion to the analysis presented by IBakos et al.l 
(2010). The HATNet photometry has already been EPD- 
and TEA-corrected before the global modeling, so we 
only considered corrections for systematics in the follow- 
up light curves. We chose the "ELTG" method, i.e., 
EPD was performed in "local" mode with EPD coeffi- 
cients defined for each night, and TFA was performed in 
"global" mode using the same set of stars and TFA coef- 
ficients for all nights. The five EPD parameters were the 
hour angle (representing a monotonic trend that changes 
linearly over time) , the square of the hour angle (reflect- 
ing elevation), and the stellar profile parameters (equiv- 
alent to FWHM, elongation, and position angle of the 
image). The functional forms of the above parameters 
contained six coefficients, including the auxiliary out-of- 
transit magnitude of the individual events. For each 
system the EPD parameters were independent for all 
nights, implying 12, and 24 additional coefficients in the 
global fit for HAT-P-18b and HAT-P-19b respectively. 
For the global TFA analysis we chose 20 template stars 
for each system that had good quality measurements for 
all nights and on all frames, implying an additional 20 
parameters in the fit for each system. In both cases the 
total number of fitted parameters (43, and 49 for HAT- 
P-18b and HAT-P-19b respectively) was much smaller 
than the number of data points (422, and 438, counting 
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only RV measurements and follow-up photometry mea- 
surements) . 

The joint fit was performed as described in lBakos et al.l 
(2010). We minimized \ 2 m the space of parameters by 
using a hybrid algorithm, combinin g the down hill sim- 
plex method (AMOEBA; see iPress et all 119921 ) with a 
classical linear least squares algorithm. Uncertainties 
for the parameters were derived applyi ng the M arkov 
Chain Monte-Carlo method (M CMC, see Ford 2006) us- 
ing "Hyperplane-CLLS" chains (]Bakos et al.ll2010t l. This 
provided the full a posteriori probability distributions of 
all adjusted variables. The a priori distributions of the 
parameters for these chains were chosen to be Gaussian, 
with eigenvalues and eigenvectors derived from the Fisher 
covariance matrix for the best-fit solution. The Fisher 
covariance matrix was calc ulated ana lytically using the 
partial derivatives given by iPail (|2009| ) . 

Following this procedure we obtained the a posteriori 
distributions for all fitted variables, and other quantities 
of interest such as a/R+. As described in Section [3TT1 
a/R± was used together with stellar evolution models to 
infer a theoretical value for logg* that is significantly 
more accurate than the spectroscopic value. The im- 
proved estimate was in turn applied to a second iter- 
ation of the SME analysis, as explained previously, in 
order to obtain better estimates of T e ff* and [Fe/H] . The 
global modeling was then repeated with updated limb- 
darkening coefficients based on those new spectroscopic 
determinations. The resulting geometric parameters per- 
taining to the light curves and velocity curves for each 
system are listed in Table M 

Included in each table is the RV "jitter". This is a 
component of noise that we added in quadrature to the 
internal errors for the RVs in order to achieve x 2 /dof = 1 
from the RV data for the global fit. It is unclear to what 
extent this excess noise is intrinsic to the star, and to 
what extent it is due to instrumental effects which have 
not been accounted for in the internal error estimates. 

The planetary parameters and their uncertainties can 
be derived by combining the a posteriori distributions 
for the stellar, light curve, and RV parameters. In this 
way we find masses and radii for each planet. These and 
other planetary parameters are listed at the bottom of 
Table IE We find: 

• HAT-P-18b - the planet has mass M p =0.197 ± 
0.013 Mj, radius R p =0.995 ± 0.052 Rj, and mean 
density p p =0.25 ± 0.04 gcm~ 3 . 

• EAT-P-1% - the planet has mass M p =0.292 ± 
0.018 Mj, radius R p =1.132 ± 0.072 Rj, and mean 



density p p =0.25 ± 0.04 gem 3 . 

Both planets have an eccentricity consistent with zero 
(e =0.084 ± 0.048 for HAT-P-18b, and e =0.067 ± 0.042 
for HAT-P-19b). As mentioned above, for HAT-P-19, 
the RV residuals from a single-Keplerian orbital fit ex- 
hibit a linear trend in time. We therefore included an 
acceleration term to account for this trend. We find 
7 =0.439±0.048ms _1 d _1 . In SectionHwe consider the 
implication of additional bodies (stellar or planetary) in 
the HAT-P-19 system. 

4. DISCUSSION 



Figure [TO] compares HAT-P-18b and HAT-P-19b to 
other known TEPs on a mass-radius diagram. We dis- 
cuss the properties of each planet in turn. 

4.1. HAT-P-18b 

From the iFortnev et al.l (]2007l ) planetary models, the 
expected radius for a coreless 0.20 ± 0.01 Mj planet 
orbiting a 4.5 Gyr star with a Solar-equivalent semi- 
major axis of 0.1073 ± 0.0071 AU is - 1.02 Rj, which 
is consistent with the measured radius for HAT-P-18b of 
1.00 ± 0.05 .Rj. The preferred age for HAT-P-18 from the 
YY isochrones (12.4+g' 4 Gyr), is somewhat older than 
4.5 Gyr, in which case the expected planetary radius 
would be even smaller. If a slight core of 10 M e is as- 
sumed, the expected radius of 0.91 Rj is below the mea- 
sured radius. We conclude therefore that HAT-P-18b is 
a predominately hydrogen- helium gas giant planet, and 
does not possess a significant heavy element core. 

HAT-P-18b is perhaps most similar in properties to 
the slightly higher density planet HA T-P-12b (M = 
0.211 ± 0.012 Mj, R = 0.959±g;gl?i2j; IHartman et all 
120091 ). Both planets orbit K dwarfs (HAT-P-18 has 
mass AU = 0.77 ± 0.03 M Q , and HAT-P-12 has mass 
M* = 0.73 ± 0.02 M Q ). However HAT-P-18 appears 
to be older than HAT-P-12, having an isochrone age of 
12.4±6.t Gyr compared with 2.5 ± 2.0 Gyr for HAT-P-12. 
HAT-P-18 is also more metal rich ([Fe/H]=+0.10±0.08) 
than HAT-P-12 ([Fe/H]=-0.36 ± 0.04). 

4.2. HAT-P-19b 

Like HAT-P-18b, HAT-P-19b also does not appear 
to possess a s i gnfica nt heavy element core. From the 
IFortnev et al.l (|2007| ) planetary models, the expected 
radius for a coreless 0.29 ± 0.02 Mj planet orbiting a 
4.5 Gyr star with a Solar-equivalent semi-major axis of 
0.0763 ± 0.0065 AU is - 1.02 Rj, which is lower than the 
measured radius for HAT-P-19b of 1.13 ± 0.07 Rj. If an 
age of 1.0 Gyr is assumed, the expected planet radius in- 
creases to 1.06 Rj, but is still slightly lower than, though 
consistent with, the measured radius. 

Like HAT-P-18b, HAT-P-19b is also very similar in 
mass/radius to another TEP, in this ca se WASP-21b 
(M = 0.30 ± 0.01 Mj, R= 1.07 ± 0.05 -Rj; IBouchv et"afl 
2010). WASP-21b orbits a somewhat hotter star than 
HAT-P-19b (WASP-21 has M = l.Olig;^Af , while 
HAT-P-19 has M = 0.84 ± 0.04M Q ). HAT-P-19 is also 
more metal rich than WASP-21 ([Fe/H]=+0.23±0.08 for 
HAT-P-19, while [Fe/H]=-0.4± 0.1 for WASP-21). 

As noted in Section [131 the RV residuals of HAT-P- 
19 show a linear trend in time, which is evidence for 
a third body in the system. The evidence for this trend 
comes entirely from the Keck/HIRES observations which 
span 144 days. Although the Subaru/HDS observations 
predate the Keck/HIRES observations, the uncertain RV 
zero-point difference between these two datasets prevents 
us from comparing the Subaru/HDS observations with 
the Keck/HIRES observations to extend the baseline for 
measuring the linear variation in the RV residuals. The 
Subaru/HDS observations span only 3 days, so the trend 
is not evident in this dataset. Following I Winn et al.l 
(|2010f ). we set 7 ~ GM C sin i c /a% to give an order-of- 
magnitude constraint on the third body, assuming the 



HAT-P-18b, and HAT-P-19b 



15 



TABLE 9 
Orbital and planetary parameters for HAT-P-18b and HAT-P-19B a 

Parameter HAT-P-18b HAT-P-19b 

Light curve parameters 

P (days) 5.508023 ± 0.000006 4.008778 ± 0.000006 

T c (BJD) b 2454715.02174 ± 0.00020 2455091.53417 ± 0.00034 

T14 (days) b 0.1131 ± 0.0009 0.1182 ± 0.0014 

T12 = T 34 (days) b 0.0150 ± 0.0008 0.0172 ± 0.0014 

a/H,, 16.04 ±0.75 12.24 ±0.67 

C/R* 20.36 ±0.08 19.76 ±0.12 

Rp/R* 0.1365 ±0.0015 0.1418 ± 0.0020 

b 2 105+ 0040 163+ 0055 

U - iUO -0.040 U - iD,;i -0.057 

b = a cos i/R+ 0.324l°'^5 0.404±°°f * 

* (deg) 88.8 ± 0.3 88.2 ± 6.4 

Limb-darkening coefficients c 

a, i (linear term, i filter) 0.4372 0.4135 

b, i (quadratic term) 0.2276 0.2459 

a,g 0.8477 0.8016 

b,g -0.0060 0.0368 

RV parameters 

Kims- 1 ) 27.1 ±1.6 42.0 ±2.1 

ifms-M" 1 ) ... 0.439 ±0.048 

fc RV d -0.035 ±0.038 -0.009 ±0.029 

h KY d 0.063 ±0.062 -0.058 ± 0.054 

e 0.084 ± 0.048 0.067 ± 0.042 

lo (deg) 120 ± 56 256 ± 77 

RV jitter (ms -1 ) 5.0 6.7 

Secondary eclipse parameters 

T s (BJD) 2454717.65 ± 0.13 2455093.515 ± 0.074 

T s 14 0.127 ± 0.014 0.107 ± 0.010 

T M 2 0.0173 ±0.0029 0.0149 ± 0.0020 

Planetary parameters 

M p (M T ) 0.197 ± 0.013 0.292 ± 0.018 

Rp (fij) 0.995 ±0.052 1.132 ±0.072 

C(Mp,R p ) ° 0.19 0.35 

p p (g cm" 3 ) 0.25 ± 0.04 0.25 ± 0.04 

loggj, (cgs) 2.69 ±0.05 2.75 ± 0.05 

a (AU) 0.0559 ± 0.0007 0.0466 ± 0.0008 

T cq (K) 852 ± 28 1010 ± 42 

f 0.029 ±0.002 0.028 ±0.002 

(F) (W 8 cvgs- 1 cm- 2 ) ^ 1.19 ±0.16 2.35 ± 0.41 

a We list the median value of each parameter from its MCMC a posteriori distribution. We 
also provide the upper and lower la error-bars about the median for each parameter. 

T c : Reference epoch of mid transit that minimizes the correlation with the orbital period. 
It corresponds to N tr — 24. T14: total transit duration, time between first to last contact; 
T12 — T34: ingress/egress time, time between first and second, or third and fourth contact. BJD 
is calculated from UTC. 

c Values for a quadratic law, adopted from the tabulations by Claret ( 2004) according to the 
spectroscopic (SME) parameters listed in Table [8] 

The Lagrangian orbital parameters derived from the global modeling, and primarily determined 
by the RV data. 

c Correlation coefficient between the planetary mass M p and radius R p . 

f The Safronov num ber is given by & = i (Vosc/V orb ) 2 = (a/ R P )(M P /M„) (see 
Hansen & Barman 2007). 
s Incoming flux per unit surface area, averaged over the orbit. 

orbit is circular. This gives riastron, and time of conjunction. The object could also 

be a low-mass star with M > 90 Mj if it has a c > 10 AU. 
' M c smi c \ f a c \- 2 

M I VI AU J ~ 4-3. Core Mass-Metallicity Correlation 

As noted in the introduction, the previously known 
The time span of the RV measurements, and the lack Saturn-mass planets exhibited a suggestive correlation 
of evidence for jerk (7) in the RV residuals, lets us between core mass (or density) and host star metallic- 
put a rough limit on the third body's orbital period of ity. The two low density planets HAT-P-12b and WASP- 
P c > 2 x 144 = 288 d, or a c > 0.8 AU. This gives a rough 21b are consistent with having no core, and orbit sub- 
lower limit on the mass of the third body of M c > 0.6 Mj, solar metallicity stars. While the higher density plan- 
though this depends on the eccentricity, argument of pe- ets Kepler-9b, Kepler-9c, CoRoT-8b, WASP-29b, and 
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Fig. 10. — Mass-radius diagram of known TEPs (small filled squares). HAT-P-18b and HAT-P-19b are shown as large filled squares. 
Overlaid arc Fortncy ct al. (2007) planetary isochrones interpolated to the solar equivalent semi-major axis of HAT-P-18b for ages of 1.0 Gyr 
(upper, solid lines) and 4 Gyr (lower dashed-dotted lines) and core masses of and 10 Mq (upper and lower lines respectively), as well as 
isodensity lines for 0.4, 0.7, 1.0, 1.33, 5.5 and 11.9gcm~ 3 (dashed lines). Solar system planets are shown with open triangles. 

HD 149026b are consistent with having substantial cores, 
and orbit super-solar metallicity stars. The apparent cor- 
relation between planet core mass and host star m etal- 
licity was prev i ously noted by iGuillot et al.l (|2006l ) and 
iBurrows et al.l (|2007l ) for all TEPs known at the time 
(nine and fourteen respectively). Many of the planets 
with M > 0.4 Mj have radii that are larger than can 
be accommodated by theoretical models, so it is un- 
clear whether the inferred core masses are physically 
meaningful for these planets. Nonetheless, for pl anets 
in the mass-range 0.4-0.7 Mj lEnoch et all ()2010[ ) find 
that planet radius is inversely proportional to host star 
metallicity, which is what would be expected if the heavy 
element content of these planets (or core mass) is propor- 
tional to host star metallicity. Figure fTTI s hows the rela- 
tion b etween core mass inferred from the iFortnev et al.l 
(2007) models and stellar metallicity for planets with 
0.15 < M < 0.4 Mj, and 0.4 < M < 0.6 Mj. HAT- 
P-18b and HAT-P-19b do not follow the correlation that 
was previously seen for the other Saturn-mass planets. 
However, since the sample size of known Saturn-mass 
TEPs is still quite small, further discoveries are needed 
to illuminate the properties of planets in this mass-range. 
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Fig. 11. — Planet core mass vs. host star metallicity for planets with 0.15 < M < Q.4Mj (left) and 0.4 < M < 0.6Mj (right). The 
core mass for each planet is determined by linear interpolation within the Fortncy ct al. (2007) planet model tables for the estimated age, 
mass, and solar-equivalent semimajor axis of the planet. We adopt an age of 4.0 Gyr or 0.3 Gyr for systems with an estimated age greater 
or less than these limits, and we adopt a solar-equivalent semimajor axis of 9. 5 AU for Saturn. Planets with a negative inferred core mass 
have radii that are too large to be accommodated by the Fortncy ct al. (2007) models. In this case the core mass is linearly extrapolated 
from the models, and provides a measure for the degree to which the observed radius disagrees with the models. The location of Saturn is 
indicated by the "S" in the left plot. HAT-P-18b and HAT-P-19b do not follow the previously suggestive correlation between core mass 
and host star metallicity for Saturn-mass planets. No correlation is apparent for planets with M > 0.6 Mj. 
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